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Gravitational waves from neutron star binary inspirals contain information about the as yet un¬ 
known equation of state of supranuclear matter. In the absence of definitive experimental evidence 
that determines the correct equation of state, a number of diverse models that give the pressure 
inside a neutron star as function of its density have been constructed by nuclear physicists. These 
models differ not only in the approximations and techniques they employ to solve the many-body 
Schrodinger equation, but also in the internal neutron star composition they assume. We study 
whether gravitational wave observations of neutron star binaries in quasicircular inspirals up to 
contact will allow us to distinguish between equations of state of differing internal composition, 
thereby providing important information about the properties and behavior of extremely high den¬ 
sity matter. We carry out a Bayesian model selection analysis, and find that second generation 
gravitational wave detectors can heavily constrain equations of state that contain only quark mat¬ 
ter, but hybrid stars containing both normal and quark matter are typically harder to distinguish 
from normal matter stars. A gravitational wave detection with a signal-to-noise ratio of 20 and 
masses around 1.4M0 would provide indications of the existence or absence of strange quark stars, 
while a signal-to-noise ratio 30 detection could either detect or rule out strange quark stars with 
a 20 to 1 confidence. The presence of kaon condensates or hyperons in neutron star inner cores 
cannot be easily conhrmed. For example, for the equations of state studied in this paper, even a 
gravitational wave signal with a signal-to-noise ratio as high as 60 would not allow us to claim a 
detection of kaon condensates or hyperons with confidence greater than 5 to 1. On the other hand, 
if kaon condensates and hyperons do not form in neutron stars, a gravitational wave signal with 
similar signal-to-noise ratio would be able to constrain their existence with an 11 to 1 conhdence for 
high-mass systems. We, finally, find that combining multiple lower SNR detections (stacking) must 
be handled with caution since it could fail in cases where the prior information dominates over new 
information from the data. 

PACS numbers: 04.30.-w,04.80.Nn,04.30.Tv 


I. INTRODUCTION 

The next generation of gravitational wave (GW) de¬ 
tectors is scheduled to start collecting data in the 
next months. Initially consisting of Advanced LIGO 
(aLIGO) [1] and Advanced Virgo (AdV) [2], the network 
of detectors will be expanded in the next few years to 
include KAGRA [3] and possibly LIGO-India [3]. Al¬ 
though the main objective of the advanced detectors is 
to provide the first direct detection of GWs, the scientihc 
community has started shifting focus to an even more in¬ 
teresting question: Once we have secured the first GW 
detections, what can we infer about Nature? 

One of the most promising GW sources for extract¬ 
ing physical information is the inspiral of compact bina¬ 
ries consisting of neutron stars (NSs). Not only are they 
expected to be the most abundant GW source [5], but 
also they are well understood. During the ^ 10^ sec¬ 
onds these GWs spend in the detectors’ most sensitive 
frequency range, they can be accurately tracked with the¬ 
oretical templates. Modeling the detectors’ output with 
these templates cannot only lift the extremely weak sig¬ 
nal out of the noise (detection), but also provide us with 
estimates of the parameters of the GW (parameter esti¬ 
mation) . 

One set of parameters that affects the templates is 


linked to finite-size effects that NSs experience because 
they are extended bodies with structure. When objects 
with a finite size are subjected to the tidal field of another 
object their multipole moments are affected in a way that 
depends on the equation of state (EoS) - for barotropic 
fluids, a relation between pressure and density- of their 
matter. The densities encountered in NS interiors are ex¬ 
tremely high; in the inner core they even exceed nuclear 
densities. In this high density regime, laboratory exper¬ 
iments and observations have still to provide a definitive 
answer on the correct EoS. We, therefore, study whether 
GWs can be used to answer the following question: Given 
the detection of the NSNS quasicircular inspiral, can we 
use finite-size effects to learn about the EoS of the ex¬ 
tremely dense NS interior EHIS]? 

To leading order, finite-size effects cause the 
quadrupole moment tensor of a star Qij to be affected 
by the tidal field tensor of its companion 5^- through 
Qij = —X£ij. The constant of proportionality A is called 
the tidal deformability and it is a function of the mass 
and the EoS. This tidal interaction causes NSs to be dis¬ 
torted during the inspiral phase and torn apart before 
merger M- By the ‘inspiral phase’ we define the evolu¬ 
tion of the binary up to an orbital separation of six times 
the total mass or contact, whichever comes first. The 
plunge and final collision after this orbital separation is 
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called the ‘merger phase’. Since this merger is expected 
to happen at high frequencies of O(10^)Hz, where the 
detector noise is likely to dominate, we here focus on the 
better modeled inspiral parlQ The relatively small veloc¬ 
ities (never exceeding 0.3 times the speed of light) of this 
inspiral phase make it ideal for a post-Newtonian (PN) 
descriptior0 in most of the frequency range considered, 
where all quantities are expanded in powers of v /c [S]. 
The PN waveform, in fact, becomes less accurate near 
the merger [IlHIS]. For nonspinning NS binaries, a more 
accurate waveform is available using the effective-one- 
body (FOB) approach [20] • In this paper, we consider 
precessing NS binaries as we will explain in more detail 
below. Since the FOB waveform for precessing NS bi¬ 
naries is currently unavailable, we use the precessing PN 
waveform. 

In the PN framework, the first finite-size effect enters 
the waveform at 2PN order in the GW phase through 
spin corrections to the quadrupole moment of the ob¬ 
jects Qi ^2 IU, and here we also include the IPN cor¬ 
rection to this [22] ■ Then, at 5PN order and above, the 
tidal deformabilities A 1.2 enter the phase directly [TTI - 
[MldiKIllEa]. All these parameters are FoS-dependent, 
however, two of the authors showed that their interre¬ 
lation is approximately FoS-independent [21] US] . A lot 
of work has been put into understanding and extending 
this result |26H43j . in particular to include higher-order 
multipole moments and tidal deformabilities [sniiisHiz]- 
We can therefore use the Love-Q relation to choose Q 1.2 
in favor of Ai ^2 or vice versa from the GW phase. 

The problem of the detectability of finite size effects 
with gravitational waves has gathered a lot of attention 
in recent years. Initial studies, based on quantifying the 
differences between waveforms with different FoSs or on 
a Fisher information matrix analysis [121 da dSl 148HS2] 
suggested that aLIGO has the potential of providing use¬ 
ful information on the NS FoS. However, due to the 
expected low signal-to-noise ratio (SNR) in aLIGO de¬ 
tections and due to the strong correlations between the 
different GW parameters, the applicability of a Fisher 
study, and the conclusions derived from it, is limited [551 - 
155] . 

For this reason, several Markov-Ghain Monte-Garlo 
(MGMG) [571IM] studies have recently been carried out 
in order to address FoS detectability in a more robust 
way in the context of Bayesian inference. The first such 
study was performed by Del Pozzo et al. [59], which 
showed that a few tens of detections of moderate bright¬ 
ness can be combined to provide strong constraints on the 
FoS, though this result seems to be highly dependent on 


^ A number of studies have examined the possibility of determin- 
ing the EoS from the merger phase (see na for an example). 
However, in the absence of a full and accurate template bank of 
merger waveforms, it is not clear how one would perform a full 
data analysis study. 

^ A term of relative order (v/c)'^^ is said to be a N-PN order term. 


the mass distribution of the sources [55] . Wade et al. [TS] 
studied the effects of systematic and statistical errors in 
FoS extraction, while Lackey and Wade [5T] employed 
a more realistic parametrization of the FoS and agreed 
that a few bright sources can determine the FoS of NSs. 

With the exception of [59] and [60], all previous work 
mentioned above consisted of parameter estimation stud¬ 
ies, where the tidal deformability is treated as a system 
parameter and searched over with an MGMG analysis. 
Any FoS that predicts a value of the deformability within 
the recovered uncertainty is compatible with the results 
of the MGMG. Any FoS that does not fall in the deforma¬ 
bility error bars can be ruled out. We choose a more di¬ 
rect approach here and compare the different FoS models 
directly, an approach known as model selection. In the 
latter, the FoS itself, instead of the tidal deformability 
is treated as an independent parameter of the system 
and the analysis allows the data to select which FoS is 
preferred. Our study is unique in that we use this tool 
to perform a comprehensive study on whether we can 
extract important physical information about the compo¬ 
sition of NSs, such as the existence of exotic species. 

The different FoS models proposed in the literature 
differ not only in the NS composition they assume, but 
also in the approximate schemes they employ to solve 
the many body Schrodinger equation. These approxi¬ 
mate schemes include approaches such as the variational 
method [52], Skyrme-Hartree-Fock (SHF) models [63] . 
Brueckner-Hartree-Fock (BHF) models [64], and rela¬ 
tivistic mean field (RMF) theory [551IHS]- The NS in¬ 
ternal composition might be that the FoS is constructed 
solely with normal matter (neutrons n, protons p, elec¬ 
trons e, muons p.), or it may contain kaon condensates 
(K), hyperons (H), pion condensates (tt), or quark matter 

(Q). 

Determining that 2 FoSs with the same composition 
but different approximations are distinguishable will not 
result in any new information about Nature. On the 
other hand, determining that 2 FoSs with different com¬ 
positions and different approximations are distinguish¬ 
able must be treated with caution. Can we claim that 
the difference between the FoSs we detected is due to 
their actual physical differences, or due to their distinct 
mathematical approaches? In order to avoid this obsta¬ 
cle, we compare FoSs that employ the same, or as similar 
as possible, approximations and differ only in their inter¬ 
nal composition. 

For example, consider an FoS, which contains normal 
matter, and is constructed with the variational methocE 
and another FoS, which is also constructed with the vari¬ 
ational method, but includes both normal matter and hy¬ 
perons. When we perform a comparison on these 2 FoSs, 
we are essentially comparing a hyperon FoS and a hyper- 
onless FoS. The result of this comparison can directly be 


® Refer to Appendix|^for a description of all EoSs used here along 
with their physical content and approximation schemes. 
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translated to physical information. If the hyperon EoS is 
preferred, we have detected hyperons in NS cores assum¬ 
ing one of the 2 EoSs is correct. If the hyperonless EoS 
is preferred, we have constrained hyperons in NS cores. 

At this stage one may reasonably object: Could we 
confuse a hyperon EoS constructed with one method 
with a normal matter EoS constructed with a different 
method? This concern can be alleviated by performing a 
large number of comparisons. Concluding that a hyperon 
EoS constructed with one method is preferred over a nor¬ 
mal matter EoS with the same method is not enough. 
We need to find as many pairs of EoSs that have been 
constructed with a variety of different methods and com¬ 
pare all of them. This approach (i) ensures we have not 
confused the effects of internal composition and approxi¬ 
mation schemes, and (ii) provides us with insight on how 
hyperons, or other particles, affect EoSs in general. Ef¬ 
fectively, by comparing many pairs of EoSs, we ‘average 
out’ any effect coming from how each pair of EoSs is con¬ 
structed and isolate the effect of the common difference 
between the 2 EoSs of all pairs i.e. the exotic matter. 

Another aspect in which our analysis generalizes pre¬ 
vious studies is that for the first time we include spin- 
precession both in our simulated signal and in the tem¬ 
plates. We use the fully analytic double-precession model 
of |67j . derived under the assumption that the binary 
components have small spins. This turns out to be an 
excellent approximation for NSs in the LIGO band, as 
their dimensionless spin parameter (the spin angular mo¬ 
mentum over the mass squared) is not expected to exceed 
^0.1 IHH]. The model has already been tested before in 
a data analysis context IMIIIQ], however, in those stud¬ 
ies we stopped our analysis at 400Hz to avoid finite-size 
effects and focus instead on the measurability of the 
masses and the spins. 

The inclusion of spin-precession in the templates is cru¬ 
cial. Reference [69] showed that allowing for the systems 
to process around the orbital angular momentum can 
break degeneracies between the masses and the spins, im¬ 
proving mass extraction by about an order of magnitude. 
This improved mass determination is directly translated 
to better A extraction, making EoSs easier to distinguish. 
This is because in the context of model selection, it is 
the EoS that is a GW parameter, and not A, which is de¬ 
termined through a relation of the form A(m, EoS) (see 
Sec. IIB). We should emphasize an important distinction 
here: the effect of better mass extraction due to preces¬ 
sion has nothing to do with the actual spin magnitude 
of the signal we study. What is important, however, is 
allowing for the template we recover the signal with to 
model precessional effects. 

The main results of our analysis are summarized below. 

We find that advanced detectors will he able to place 
strong constraints on the existence of quark stars com¬ 
prised solely of quark matter. Given the EoSs available 
today, a NSNS binary with masses in the (1.2, 1 . 5 )M 0 
range with SNR = 30 can effectively rule out strange 
quark stars, or make a positive detection of them. We, 


furthermore, argue that even for the plausible EoSs 
constructed in the future, there exists some mass € 
(1 Mq,1.8Mq) where quark stars are distinguishable 
from normal matter NSs. 

The prospects of detecting or ruling out hybrid NSs 
including both normal and quark matter are worse. If 
the strong interactions between quarks are close to those 
predicted by a perturbative analysis |72j and the transi¬ 
tion between nuclear and quark matter phases happens 
around twice the nuclear saturation density, the detection 
of a (1.4, 1 . 35 )M 0 NS binary with SNR 30 — 40 could pro¬ 
vide significant evidence of whether quarks form in NS 
interiors. However, if the strong interactions are weaker, 
aLIGO will not be able to reach confident conclusions. 

It is unlikely that aLIGO will be able to claim a de¬ 
tection of hyperons or kaons, since that would require 
high masses and SNR > 60. The detection of hyperon 
or kaon condensates in NS interiors requires high mass 
stars, since it is only at these high masses that you en¬ 
counter densities large enough for these condensates to 
form. This poses a significant problem; the importance of 
finite size effects is reduced with increasing mass since A 
decreases with increasing mass. Moreover, most NSs are 
expected to have masses around ~ 1 . 4 M 0 , rarely reach¬ 
ing the 2 M 0 required for hyperons and kaons detection. 
We therefore conclude that it is unlikely that aLIGO will 
be able to positively identify hyperons or kaons in NSs. 
On the other hand, if hyperons and kaons are not formed 
in NS interiors, aLIGO could place constraints on their 
existence. 


Our analysis suggests that aLIGO can distinguish be¬ 
tween models that differ at low central densities, like nor¬ 
mal matter EoSs and EoSs containing quark matter. In 
order to probe the high density regime we need SNRs 
higher than what aLIGO is likely to achieve. This is due 
to the fact that high mass systems (i) present smaller 
finite size effects, and (ii) have masses close to the maxi¬ 
mum mass allowed, which causes some interesting effects 
related to the prior boundary (see Sec. IID). 

Among the various noise configurations aLIGO can 
be tuned to, the optimal for EoS determination is the 
default Zero-Detuned, High-Power (Zero-Det., High-P) 
one Tuned conhgurations include the NSNS Opti¬ 
mized (NSNS Opt.) and the High Frequency (High F.) 
ones; both perform in an inferior way when it comes to 
EoS extraction. In the case of NSNS Opt. this is due to 
its low sensitivity at frequencies above 600Hz, when it 
is exactly at these frequencies that finite size effects are 
more prominent. As far as High F. is concerned, its im¬ 
proved sensitivity is limited to a very narrow frequency 
range around 10^Hz. This fact coupled to its worsened 
sensitivity at low frequencies makes High F. unsuitable 
for EoS studies. We conclude that in the high frequency 
regime we are interested in, the Zero-Det., High-P. noise 
curve has the overall higher sensitivity and more accu¬ 
mulated SNR. 


Stacking, i.e. combining multiple low SNR sources, 
might improve the results obtained here, but it might also 
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lead us to erroneous conclusions about the true EoS. We 
find that when the SNR is low our results can be dom¬ 
inated by prior information, rather than any new infor¬ 
mation we get from the GW data. This could lead to 
each individual low SNR binary system providing some 
confidence, but in favor of the wrong EoS. Stacking all 
these systems will inevitably lead to great confidence in 
favor of the wrong EoS. We emphasize that stacking must 
always be treated with caution. 

As far as the spin of the bodies is concerned, we find 
that the magnitude of the injected spin angular momen¬ 
tum has a negligible effect on our results, provided that 
the templates allow for spin precession. This is in agree¬ 
ment with Ref. [69) . where the order-of-magnitude im¬ 
provement in mass extraction was achieved over spin- 
aligned templates even for nonspinning systems, as long 
as spin-precessing templates were used in the recovery of 
the signal. 

The remainder of the paper is organized as follows. In 
Sec. we describe the techniques and simulations we 
use. In Sec. Ill we describe in more detail the results of 
the EoS comparison. In Sec. |IV[ we summarize our work. 
Throughout the paper we use units where G = c = I. 


II. MODEL SELECTION 

Model selection in the Bayesian framework requires an 
explicit statement of the models compared. In this sec¬ 
tion, we describe in detail the models we use, as well 
as our methodology when comparing them. We describe 
the Bayes factor (BE), a quantity that assesses which 
model is preferred by the data, and give an overview of 
our simulated signals. We conclude this section with a 
brief discussion of the power of stacking signals versus 
detecting a single loud signal. 


A. Bayesian inference 


In the context of Bayesian inference, the probability 
that a hypothesis Hi is correct given some data d is m 

[SH] 


p{Hi\d) 


p{Hi)p{d\Hi) 

p{d) 


( 1 ) 


where p{Hi\d) is the posterior belief in the hypothesis af¬ 
ter the data has been analyzed, p{Hi) is the prior belief 
based on all information we have before analyzing the 
data, and p{d) is the probability of the data, an unim¬ 
portant normalization constant in our case. The evidence 
p{d\Hi) is given by an integral over the parameters of the 
model 0 


p{d\Hi) = J d6 p{6\Hi) p{d\eHi), ( 2 ) 

where p{6\Hi) is the prior information on the model 
parameters, and p(d\6Hi) is the likelihood, where 


\np{d\6Hi) = —l/2{s — h\s — h) in Gaussian noise, 
with s the signal, h the template model and (-j •) the 
noise-weighted inner-product [5711^. 

When we have to select between 2 competing hypothe¬ 
ses, we compare their posterior beliefs through the odds 
ratio (OR), defined by 

^ p{Hi\d) ^ p{Hi)p{d\Hi) 
pmd) p{H,)p{d\H,y 

The OR is the ‘betting odds’ of Hi compared to H 2 and 
includes the both prior belief in each hypothesis and any 
new information that is extracted from the data. These 
two contributions can be separated by defining the Bayes 
factor (BE) 


p{d\Hi) 

pid\H2)’ 


(4) 


which includes only the data contribution to the OR. 
A BE > I means that Hi is supported better by the 
data, while a BE < I means that H 2 is preferred. In 
the case of uninformative priors, i.e. p{Hi) =p{H 2 ), the 
OR equals the BE. In this paper, we choose to work with 
the BE instead of the OR because we are interested in 
whether the data lends more support to some hypothesis 
over another, irrespective of our prior beliefs in them. Of 
course, once we can confidently quantify our prior belief 
in a hypothesis, we can trivially go from the BE to the 
OR. 

Working with the BE, however, has one major prob¬ 
lem: we cannot draw the same conclusions about two 
different pairs of hypotheses that have the same BE. For 
example, consider the problem of whether a given GW 
signal is better described by GR or by a modified grav¬ 
ity theory. In this case, we have a strong prior belief in 
favor of GR, given the many successes of experimental 
relativity [ZH [75] . We would therefore require very large 
BFs in favor of the modified gravity hypothesis to claim 
a detection of a deviation from GR. On the other hand, 
if we are interested in whether the signal is better de¬ 
scribed by one of two competing EoSs, we do not have 
such strong prior beliefs in favor of any of the models. 
That means that we would not need such large BFs to 
claim that we have identified the correct EoS of Nature, 
provided we remain agnostic about the two EoSs. 

For this reason, it is crucial that we explicitly define 
what we mean by “a BE that is large enough” on a 
problem by problem basis. As we will discuss later in 
this section, we work with models for which we do not 
have strong experimental prior knowledge. Therefore, we 
adopt the Jeffreys scale of interpretation of BFs m to 
define how significant a BE is. When BE < 1 it is neg¬ 
ative, for 1 < BE < 3 it is barely worth mentioning, for 
3 < BE < 10 it is strong, for 10 < BE < 100 it is very 
strong, and finally for BE > 100 it is decisive. This is 
in contrast to the analysis of Refs. [77ll80j . which dealt 
with tests of GR, where BFs around 100 were considered 
strong and not decisive, given the strong prior in favor of 
GR. 
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In order to calculate the BF between two models, a 
number of different techniques can be used: thermody¬ 
namic integration [551ISI] , nested sampling [5T] and re¬ 
verse jump, Markov chain-Monte Carlo m- Here, we 
employ the third technique which requires promoting the 
model to a parameter of the Markov chains. Then, the 
BF is given by 

time the chains spend in model 1 
time the chains spend in model 2 ’ 

with error bars calculated with the technique suggested 

in [52] ■ 


B. Models 


Our models represent GWs emitted in the late inspi¬ 
ral of NS binaries assuming GR is correct from the time 
they enter the LIGO band until the NSs come into con¬ 
tact, or their separation becomes six times the total mass, 
whichever comes first (see IIG for details). The differ¬ 
ence in the models will only be in the finite size effects 
they include. Gomparison between these models can be 
viewed as traditional parameter estimation, only now the 
model itself is an extra discrete parameter. In parameter 
estimation, the waveform of a NS binary inspiral with a 
circular orbit depends on the following parameter^ 

0PB = {mi,1712, On, 4>n, Dl,9l, 4>L,tc,4>c, Si, S2,\i,\2\, 

( 6 ) 


while in a model selection study the template depends 
on the parameters 


^MS — {Wi, m2, 4>n, Dl, 9L,4'L,tc, 4’c, Si, S2, EoS}, 

(7) 

where are the component masses, {9n,4>n) give the 
sky location of the source, is the luminosity distance, 
{9l, 4>l) give the direction of the orbital angular momen¬ 
tum, {tc, 4>c) are the time and phase of coalescence re¬ 
spectively, Si = Xim^(sin cos (/)i, sin 0 i sin (/)i, cos 0 i) are 
the spin angular momentum vectors of each binary com¬ 
ponent, with Xi ■= \Si\/m^ the dimensionless spin pa¬ 
rameter, and Xi = Xijvcvl are the dimensionless tidal 
deformabilities. The two parameter sets are equivalent, 
since knowledge of (mi,m 2 ,EoS) from 6^^ can be used 
to construct the quantities Ai(mi,EoS) and A 2 (m 2 ,EoS) 
in ^pE. 

The EoS is what defines a model. For some EoS the 
model predicts a GW with 0pE and Ai ^2 as a function 
of (mi^ 2 ,EoS) for (mi,m 2 ) < M„,„,(EoS), and no GW 
otherwise. The quantity M,„j,„(EoS) is the maximum NS 


^ In this argument, we encode the EoS dependence through 
the dimensionless tidal deformabilities Ai_ 2 ; the argument goes 
through if one uses any other quantity to parametrize tidal de¬ 
formations. 


EoS 

Composition 

AP4 [84], SV [85], SGI [85], SkI4 [86] 
DBHF(2)(A) [HI], MPa [HI], G4 [HD 
GA-FSU2.1 [90], Shen [M]|66] 

n, p, e, ^ 

SGI-YBZ6-SAA3 [9T], N1Y5KK* [HI], 
SkI4-YBZ6-SAA3 [9T] 

MPaH [HI], H4 [H2] 

n, p, e, fi, H 

SGI178 [H5] 

SV222 [HS] 

GA-FSU2.1-180 [901 

n, p, e, fj., K 

ALF4 [92], ALF5 

n, p, e, fj., TT, Q 

SQM3 [93] 

Q (u, d, s) 


TABLE I. Classification of EoSs with respect to internal com¬ 
position. The hrst cluster corresponds to EoSs with normal 
matter. The second and third clusters include hyperons and 
kaon condensates respectively. The last two rows list EoSs 
that include quark matter. 


mass that can be stably supported, given an EoS. For the 
GW template itself, we will use the small-spins, double- 
precessing waveform of EZ], which is constructed in the 
Fourier domain through a stationary-phase approxima¬ 
tion in a post-Newtonian expansion. We use an extended 
version of these templates by adding finite-size effects at 
2 and 3PN order (due to the quadrupole moment) pTl[5^ 
and at 5 PN order and higher (due to tidal deformability 
effects) [TTUISl l23| . More specifically, we include tidal 
terms that depend on the £ = 2 electric tidal deformabil¬ 
ity A up to 7.5PN order given in m- We also include the 
contribution from the i = 2 magnetic tidal deformability 
a 2 at 6 PN order and the £ = 3 electric tidal deformabil¬ 
ity A 3 at 7PN order [JS]. We further take into account 
the leading correction to the adiabatic approximation en¬ 
tering at 8 PN order that depends on the £ = 2 f-mode 
frequency /2 of a NS [HI [12]. (72, A 3 and /2 can approxi¬ 
mately be expressed in terms of A thanks to the universal 
relations found in nans]- In addition, we use the Love- 
Q relations [241125j to rewrite the quadrupole moment in 
terms of the tidal deformability. 

Which EoS models should we allow the data to select 
from? There is a great number of EoSs available in the 
literature, varying both in the type of matter they con¬ 
sider (quarks, hyperons, kaons, muons, pions, neutrons, 
protons, electrons), and in the approximation schemes 
used to construct the EoS (see Appendix 0 . We can 
classify the EoSs by the type of matter they employ; the 
subset of EoSs discussed in this paper are presented in 
Table |T| The EoSs in each of these categories differ in 
the approximation schemes used to solve the many-body 
Schrodinger equation (see Appendix 0. 

Given the great variety of EoSs, how should we carry 
out a model selection study? Gould we select a character¬ 
istic set of models, perhaps chosen by looking at how the 
A — TO relations behave for a set of EoSs? Figure shows 
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FIG. 1. Dimensionless tidal deformability as a function of 
mass for a number of EoSs with very different physical con¬ 
tents. The lack of clustering in the X — m space shows that we 
cannot perform model selection with all EoSs simultaneously. 


these relations for a few of the EoS listed above. Notice 
how the normal matter EoSs (AP4, Shen) form a band 
in the X — m space that contains both SQM3 (consists 
of pure quark matter) and H4 (contains hyperons). This 
lack of clustering with respect to the internal composition 
classification in the X — m space indicates that comparing 
arbitrary EoSs will not produce physically meaningful re¬ 
sults. For example, say we concluded that AP4 (red solid 
line) is distinguishable from H4 (black, dot-dashed line). 
Can we claim that hyperons are detectable from GW ob¬ 
servations? The answer is no; there is another normal 
matter EoS (Shen, turquoise dotted line) much closer to 
H4 than AP4. Unless we compare Shen to H4 too, we 
cannot claim detectability of hyperons. 

Another possibility would be to compare each model 
against every other model simultaneously, but (i) the 
computational cost would be prohibitive, and (ii) it is not 
clear what physical question such a comparison would 
address. Comparing models within each category that 
differ only in the approximation scheme used to calcu¬ 
late the EoS is analogous to comparing post-Newtonian 
models for the GWs emitted in the inspiral of compact 
objects. In the end, all of such models are only approxi¬ 
mations to the exact solution of Nature. Instead, we are 
interested in specific physical questions regarding the in¬ 
ternal composition of NSs. To do so, we will choose pairs 
of EoSs that are as similar as possible (using the same 
approximation schemes), but differ only in the matter de¬ 
grees of freedom. In this way, we can directly translate 
distinguishability between models into distinguishability 
between physical scenarios. 


C. Signal injections 

Whether the data can distinguish a given EoS model 
over another depends strongly on the particular sig¬ 
nal detected. The parameters that affect EoS dis¬ 
tinguishability the most are the two masses, the dis¬ 
tance to the source, which effectively controls the 
SNR, and the EoS itself. We therefore simulate dif¬ 
ferent signals with various values for these parame¬ 
ters. For the remaining parameters we select the fol¬ 
lowing injection values: (cos(/lAr) = (—0.105,3.705), 
{cosOl,4>l) = (0.801,3.216), (U,(/>c) = (1,024s,4.461), 
(cos0i,(()i) = (0.774,2.248), (cos 6 » 2 ,<(> 2 ) = (0.968,5.311), 
and (xi,X 2 ) = (0.04,0.04). All parameters have been 
randomly chosen so that they do not lead to any ‘spe¬ 
cial’ orientation of the binary (face on, edge on etc.). We 
have also performed simulations with other dimension¬ 
less spin magnitudes and found that the spin has a very 
small effect on EoS distinguishability. 

The initial misalignment of the spin and the orbital 
angular momenta means that the system will undergo 
precession [94]. Indeed, with these choices of parame¬ 
ters, the angle between the orbital angular momentum 
and the total spin angular momentum is ~ 30°. Spin 
precession is modeled through the small-spins, double¬ 
precession approach of |67|, which has been shown to be 
highly accurate for modeling NS binaries ITD). This ap¬ 
proach is valid in the inspiral phase only, since it is based 
on a post-Newtonian expansion. The signal is thus mod¬ 
eled with the same approach as the templates that define 
the EoS models. 

The approximations used to describe the orbital mo¬ 
tion in the small-spins, double-precession approach are 
only valid up to a given frequency. We thus carry out 
our analysis up to min(/isco, /c(EoS)), where /isco is 
the Keplerian frequency at r = 6 (mi -|- m 2 ) [S3 ISSj and 
fc is the contact freguency, that at which the separation 
of the two bodies is equal to the sum of their radii. Al¬ 
though Mandel et al. m showed that terminating the 
waveforms at a certain frequency can affect the results 
through the addition of artificial information, our cutoffs 
are at such high frequencies that they are not expected 
to affect the results. 

Even though the particular noise realization in the de¬ 
tector at the time the GW passes through will have an 
effect on parameter extraction, we do not inject noise 
in our analysis. Given that it is impossible to predict 
the noise instance, the best we can do is average our 
results over multiple noise realizations. However, Nis- 
sanke et al. [15] showed that such averaging is equivalent 
to zero injected noise in the signal. Sampson et al. m 
showed that a given noise realization causes the likeli¬ 
hood to shift as a whole, without significantly changing 
its shape (Fig. 4 of [75]), suggesting that the integral 
of the likelihood (the evidence) is minimally affected by 
noise fluctuations. This picture, however, does not hold 
for poorly constrained parameters, like the ones studied 
here, where the posterior extends over a large fraction of 
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the prior volume musa. 

Of course, the noise curve of the detector does play a 
very important role in the calculation of BFs, through 
the noise-weighted inner product in the likelihood m- 
We will mostly adopt the zero-detuned, high-power noise 
spectral density of the advanced detectors m, though 
we will explore other choices in Sec. |IIIE[ The specific 
form of the likelihood we use assumes that the noise 
in the detectors is stationary and Gaussian [see below 
Eq. ([^], neither of which is strictly true. Cornish and 
Littenberg, however, have shown how to model the non- 
Gaussian features |82j , and deal with the nonstationarity 
of the noise [M], leaving us with only the modeling of the 
stationary, Gaussian noise component. 

When recovering the parameters, we use a uniform 
prior in (O.1,3)M0 for the masses, a uniform prior on 
the sphere for all directions, a uniform prior in (0,1) for 
the dimensionless spin magnitudes, and a uniform prior 
in the log of the distance. All prior ranges are selected 
such that they are wide enough to not affect our results. 


D. Stacking vs high SNR 


Comparison 

EoSs 

mi(M0) m2(M0) 

kaons 

SV/SV222 

1.4 1.35 
1.95 1.9 
2.05 2 

GA-FSU2.1/GA-FSU2.1-180 

1.9 1.8 

1.99 1.95 

hyperons 

G4/H4 

1.4 1.35 

1.8 1.7 

1.95 1.9 

MPa/MPaH 

1.95 1.9 

2.15 2.1 

quarks 

ALF-GCR/GCR 

1.5 1.4 

ALF5/AP4 

1.4 1.35 

1.8 1.7 

SQM3/AP4 

1.2 1.1 


1.4 1.35 

1.8 1.7 


TABLE II. Simulated masses for the comparison of Sec. |III| 
The first column gives the particle whose existence the com¬ 
parison constrains, the second column gives the EoSs com¬ 
pared, while the third and fourth give the masses. 


In reality aLIGO will probably reach physically inter¬ 
esting conclusions by combining information from multi¬ 
ple detections, rather than by waiting for a very loud one. 
Along those lines, one would argue it makes more sense to 
stack a sufficient number of moderate SNR sources rather 
than study the BF as a function of the SNR. Our results, 
however, suggest that stacking should be performed with 
caution. 

For low SNR detections we find 2 rather counterintu¬ 
itive effects: (i) it is possible for the wrong model to be 
preferred over the correct one, and (ii) it is possible for 
the correct model to be preferred less and less as the SNR 
increases (see Appendix . The first effect is not new; 
it has already been encountered in the context of com¬ 
paring models with different dimensionality where it 
takes the form of an Occam Penalty on the more compli¬ 
cated model. The second effect is perhaps less familiar: 
why would it be that, as the signal strength increases, 
the data fails to increasingly support the correct model? 
The answer to this question and the root cause of these 
effects can be traced down to sharp cutoffs of the prior 
distribution. 

When computing BFs, we must compare models with 
different maximum masses. In fact, most of the time 
the difference is rather large, with one model allowing 
NSs up to ~ 2M0, while the competing model going up 
to ~ 2.5Mq. When the injected mass is close to the 
maximum allowed mass and the SNR is sufficiently low, 
the posterior distribution is affected by this cutoff. In 
order to understand and visualize this effect, we construct 
a simple 1—D model in Appendix [C| In the context of this 
simple model we explain both the effect on the BFs that 
favor the wrong model, and the BFs that decrease with 
increasing SNR. 


These results show that stacking many weak sources 
is not necessarily equivalent to a single bright source. 
When each new observation is informative, i.e. the likeli¬ 
hood dominates over the prior, then the data will prefer 
the correct model and BF > 1. In that case, adding the 
extracted information from this observation will push the 
analysis in the right direction, and eventually, we will re¬ 
cover the same results as from a single loud event. On 
the other hand, if the observation is not informative and 
the result is prior dominated rather than likelihood dom¬ 
inated, the wrong model might be preferred and BF < 1. 
In that case, adding this observation in the stack will 
lead the analysis in the wrong direction: a large num¬ 
ber of weak observations that favor the wrong model will 
build confidence in the wrong conclusion. 

The above results suggest that stacking must be 
treated with caution. For sufficiently high SNR events, 
the posterior will be narrow enough that it will not be af¬ 
fected by the maximum mass cutoff. In this case, one re¬ 
covers the expected result: the correct model is preferred 
and it is preferred more and more as the signal strength 
increases. However, if the signals are of lower SNRs, the 
observations may not be informative, and then, the final 
cumulative result may be largely influenced by the prior 
and not by the new information contained in the signals. 


III. COMPARING EQUATIONS OF STATE 


Our goal is to study whether GW detections of in- 
spiraling NSs can be used to learn about their interior 
composition and in particular, whether they contain kaon 
condensates (Sec. Ill A), hyperons (Sec. Ill B), and quark 



























matter (Sec. IIIC). To do so, we need to isolate their re¬ 
spective effects in the EoSs. We accomplish this by com¬ 
paring pairs of EoS models that are as similar as possible, 
but differ only in the inclusion of one of these particles. 
Appendix [^provides a comprehensive classification of all 
the EoSs we use in this section. 


A. Kaon condensates 


We first address the question of whether NSs with kaon 
condensates in their inner cores leave an observable sig¬ 
nature on inspiral GWs. To do so, we choose 3 pairs 
of EoSs constructed with (i) the Skyrme-Hartree-Fock 
(SHF) scheme, (ii) relativistic mean field (RMF) theory, 
and (hi) the SHF scheme including three-nucleon interac¬ 
tions (TNI)[^ Each pair consists of one EoS with a kaon 
condensate and one without. 

Kaon condensates could emerge in stars at high central 
densities and, therefore, kaon models differ from kaon- 
less ones only for NSs with sufficiently high masses. This 
makes the extraction of physical information from these 
systems more difficult than that for low-mass systems for 
3 reasons. From a data analysis point of view, systems 
with masses close to the maximum allowed mass will suf¬ 
fer from the edge effects described in Sec. HD[ making it 
more difficult to get likelihood-dominated results. From a 
physical point of view, NSs with high masses have smaller 
values of A, as seen in Fig.[^ making finite-size effects less 
relevant in the GW phase. From an astrophysical point 
of view, NSs with masses around 2Mq are expected be 
rare. Table [TTl the values of the masses we select for each 
EoS comparison. 

As we will show below, we find that if kaon condensates 
indeed form in the inner core of NSs, they will be hard 
to detect with the signals expected from aLIGO. On the 
other hand, if condensates do not form in NSs, then loud 
aLIGO signals may be able to establish this with a high 
mass detection. 



FIG. 2. m-R relation for the EoS pairs that test kaons: 
The black lines correspond to SV (solid) and SV222 (dotted), 
which are constructed through SHF models. The red lines are 
for GA-FSU2.1 (solid) and GA-FSU2.1-180 (dotted), which 
are constructed through RMF theory. The presence of kaons 
causes the kaon model of each pair to differ from the normal 
matter model for high masses. The pairs of similar blue sym¬ 
bols indicate the component masses we use in our analyses. 


3. SHF models + TNI |85j : 

(a) kaonless: SGI. 

(b) with kaons: SGI178. 

The results of comparison (3) are very similar to the re¬ 
sults of comparison (2), so we will only present pairs (1) 
and (2). Figure]^ shows the mass-radius relation for the 
EoSs we present. This figure suggests that kaons affect 
the properties of only the most massive NSs. The pairs 
of circles, triangles, and squares indicate the values of the 
masses in the injected signals of our analysis. 


1. Models 

The 3 pairs of models we compare are constructed with 
different approximations. 

1. SHF models |85| : 

(a) kaonless: SV. 

(b) with kaons: SV222. 

2. RMF theory [50] : 

(a) kaonless: GA-FSU2.1. 

(b) with kaons GA-FSU2.1-180. 


® See Appendix |A| for a short description of each method. 


2. Bayes factor 

Figure shows the BF in favor of the kaon model of 
each pair (left panel) and the kaonless model of each pair 
(right panel) as a function of the SNR of the signal. The 
different lines correspond to the different injected masses 
indicated by the blue symbols in Fig. We always plot 
the BF in favor of the correct (injected) model. Any BF 
> 1 in the left panel means that we correctly detected 
the presence of kaons in the interior of the NSs, while 
any BF > 1 in the right panel means that we correctly 
concluded that there are no kaons in their interiors. 

As first suggested by the mass-radius plot. Fig. [^con¬ 
firms that in order to detect the presence of kaons on 
NS EoS we need a high mass system. The lowest mass 
system (mi = 1.4M0,m2 = 1.35M0) gives BFs that are 
consistent with 1, in agreement with the m—R relation. 
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FIG. 3. BF in favor of the kaon models (left) and the kaonless models (right) as a function of the SNR for different injected 
masses given in the brackets. BFs are always quoted in favor of the correct (injected) model. For example, the black dotted- 
dashed curve labeled “SV222 (1.4,1.35)” in the left panel means that the injected model is the SV222 EoS with the NS mass 
(1.4,1.35)M0, and Hi = SV222 and H 2 = SV in Eq. BF > 1 in the left panel means that we correctly detected the 
presence of kaons inside the NSs, while that in the right panel means that we correctly ruled out such a presence of kaons. The 
kaonless models give higher BFs than the kaon models, making it easier to establish that NSs do not have kaon condensates 
than the opposite. BF < 1 in the low SNR regime for high-mass systems in the left panel is due to the fact that the injected 
masses are close to the maximum mass allowed in each model (see the main text for more detail). 


For large SNRs, as we increase the injected masses we re¬ 
cover BFs that correctly favor the kaon model (left panel) 
and correctly favor the kaonless models (right panel). 

At low SNRs, however, we encounter BFs that incor¬ 
rectly disfavor the kaon models on the left panel of Fig. 
This is because the injected masses are very close to the 
maximum mass allowed by the correct model. When an 
injected parameter is closer to the edge of the prior than 
the width of the posterior, the posterior has to be cut off 
(see Sec. IID and Appendices [b] and [C|) . In our case, the 
injected masses are close enough to the maximum mass 
of the kaon models that these models “lose” some chain 
points because the models simply cannot produce such 
high mass systems. What the kaon models are in fact 
trying to do when the mass is above their maximum al¬ 
lowed mass is match the GW signal with pure Gaussian 
noise. How this affects the BF is clear from Eq. 0. Re¬ 
call that the BF in favor of SV222, for example, is equal 
to the number of times the chains visited this model, di¬ 
vided by the number of times they visited the competing 
model SV. When SV222 has an abrupt cutoff close to 
the injected masses, a lot of the chain points will be dis¬ 
favored because the model fails to produce a signal for 
these values of the masses. This leads to the chains visit¬ 
ing the SV222 model less, and in the end a BF that does 
not favor the correct model. 

The discussion of Sec. lIIDI indicates that when we have 
BFs favoring the wrong model, we are in a regime where 
it is mostly the prior that dominates our results. Obvi¬ 
ously, the results in this regime do not offer new physical 
information and cannot be used to claim that we con¬ 


strained the EoS. For example, imagine we detected a 
(1.99, 1.95)Mq system at SNR 30 and we want to claim 
something about the presence or absence of kaons in the 
system. If we choose to test whether GA-FSU2.1 or 
GA-FSU2.1-180 fit the data better, we will recover BFs 
in favor of the kaonless model GA-FSU2.1 regardless of 
whether kaons are actually present or not. The data com¬ 
ing from such a system can clearly not be trusted to give 
the correct result. 

For the case studied here, the kaon model comparison 
starts to become likelihood-dominated (BF > 1) when 
the SNR > 40 in the SV/SV222 case, and when the SNR 
is somewhere above 60 in the GA-FSU2.1/GA-FSU2.1- 
180 case. A detection of a (2.05, 2)Mq system at SNR = 
40 and above will lead to a strong BF in the Jeffreys scale 
in favor of the kaonless model for a kaonless injection. On 
the other hand, for a kaon injection, the SNR needs to 
reach 60 and above in order to obtain strong BFs. For 
lower mass systems, the BFs are lower, and barely worth 
mentioning in the Jeffreys scale. 


B. Hyperons 

Moving on to the study of the detectability of hyperons 
in the inner cores of NSs, we select 5 EoSs constructed 
with (i) RME theory, (ii) a nonrelativistic Brueckner- 
Hartree-Fock (BHF) approach, (iii) a relativistic BHF 
approach, (iv) and (v) a SHF approach. Each pair con¬ 
sists of one EoS with a hyperon in the inner core and one 
EoS without. 
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FIG. 4. m-R relation for the EoS pairs that test hyperons: 
The black lines correspond to G4 (solid) and H4 (dotted), 
which are constructed through RMF theory. The red lines 
are for Mpa (solid) and MPaH (dotted), which are constructed 
through the nonrelativistic BHF method. The presence of hy¬ 
perons causes the hyperon model of each pair to differ from 
the normal matter model for high masses. The pairs of sim¬ 
ilar symbols indicate the component masses we use in our 
analyses. 


Hyperons, much like kaons, affect the EoS of only the 
most massive NSs, since they form at the very highest 
central densities. For that reason, the results of this sec¬ 
tion are very similar to the previous one on kaons: we 
find that aLIGO can constrain the existence of hyper¬ 
ons in NSs, but detecting them will be much harder (see 
Table |!!|. 

1. Models 

The models we use to determine hyperon detectability 
are the following. 

1. RMF theory [8^ : 

(a) hyperonless: G4. 

(b) with hyperons: H4. 

2. Nonrelativistic BHF models |88j : 

(a) hyperonless: MPa, 

(b) with hyperons: MPaH. 

3. relativistic BHF models m- 

(a) hyperonless: DBHF(^)(A). 

(b) with hyperons: N1Y5KK*. 

4. SHF model: 

(a) hyperonless: SGI [55] . 


(b) with hyperons: SGI-YBZ6-SAA3 [HI]. 

5. SHF models: 

(a) hyperonless: SkI4 [86] . 

(b) with hyperons: SkI4-YBZ6-SAA3 [91] . 

The results of comparisons (3), (4) and (5) are very sim¬ 
ilar to the results of the MPa and MPaH comparison, so 
we will not present them here. Figure]^ shows the mass- 
radius relation for the EoSs for which we will present 
comparisons. Glearly, hyperons affect the EoS of only 
the most massive NSs, exactly like kaons. The pairs of 
symbols indicate the masses we inject. 


2. Bayes factors 

The results of our analysis are presented in Fig. [^ 
which shows the BF in favor of the hyperon models (left) 
and the hyperonless models (right) as a function of the 
SNR of the injected signal. The message of this plot is 
clear if we take into account the reasoning presented in 
the previous section. In the G4-H4 comparison for the 
hyperonless comparison at low SNR (below 40), we en¬ 
counter the effect of the wrong model being preferred over 
the correct one. On the other hand, in the MPa-MPaH 
comparison, this continues to hold until the SNR = 60. 
This indicates that it is the prior (and more specifically 
the maximum mass of the hyperon model) that domi¬ 
nates our results and not the likelihood. When the SNR 
exceeds 40, we start extracting interesting information 
from the comparison. Assuming hyperons do not form in 
NS cores, a (1.95, 1.9)Mq detection with SNR = 40 and 
60 will give a strong and very strong indication in the 
Jeffreys classihcation of the nonexistence of hyperons re¬ 
spectively. On the other hand, if hyperons do form in NS 
cores, then a signal with SNR = 60 would only provide 
strong evidence, as the BFs do not exceed 7. 


C. Quark matter 

Unlike kaon condensates and hyperons that can only 
exist in combination with normal matter, quark matter 
can form both with and without normal matter. In the 
first case we have hybrid NSs with EoSs of the ALF [55] 
type that have quarks formed after a certain transition 
density. The second case results in quark stars [53] com¬ 
prised solely of quark matter. 

The EoSs of pure quark stars differ from normal matter 
EoSs even at low densities, making it possible for aLIGO 
to detect or rule them out. On the other hand, hybrid 
normal/quark matter EoSs are constructed by stitching 
the nuclear matter EoSs in the low density regime to 
quark matter ones in the high density regime, with ap¬ 
propriate phase transitions in between. Therefore, they 
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FIG. 5. BF in favor of the hyperon models (Left) and the hyperonless models (Right) as a function of the SNR for different 
injected masses. BFs are always quoted in favor of the correct (injected) model. We conclude that it is easier to constrain 
hyperons than detect them. A similar result was reached in the case of kaons as shown in Fig. 


reduce to normal matter EoSs at low densities. The tran¬ 
sition density ric and the strength of the strong interac¬ 
tions c determine how much the hybrid EoS differs from 
the normal matter EoS it is stitched to. The constant 
c is predicted to be c ^ 0.37 m by a perturbative cal¬ 
culation. When c is close to this value and the transi¬ 
tion from a nuclear matter to quark matter happens at 
around twice the saturation density, we find that hybrid 
EoSs might be distinguishable from normal matter EoSs 
for signals with SNR ^ 30 — 40. 

1. Hybrid EoSs 

In order to study hybrid normal/quark matter stars 
we select EoSs of the ALE family [52] , the normal matter 
part of which are based on AP4 [84]. All EoSs include 
ordinary matter, while the ALE EoSs also include pions 
and quarks in the inner cor^ The left panel of Eig. 
shows the mass-radius relation for ALF5 and AP4. ALF5 
has (nc,c) = (2no,0.4), where no = 0.16 fm“^ is the 
nuclear saturation density. 

We also constructed new hybrid star EoSs (GCR-ALF) 
by stitching the nuclear matter GCR EoS constructed 
in [100] to the same quark matter EoS as the ALE 
family. The m-R relation for a GGR-ALF EoS with 
(ric, c) = (2no,0.35) is shown again in Fig. together 
with the corresponding nuclear matter GCR EoS with 
the symmetry energy of Egy^ = 33.8MeV. We have also 


® In principle, we should treat pions as a separate particle, like 
kaons and hyperons. However, we are not aware of any EoS 
model that includes only pion condensates and predicts a maxi¬ 
mum NS mass above 2Mq. 


chosen different values for uq , c and Egym and found that 
the difference between the nuclear matter and hybrid 
EoSs are typically even smaller than the one in Fig. 

The right panel of Fig. [^ shows the BFs we recover for 
the ALF5/AP4 (black) and the ALF-GCR/GCR com¬ 
parison (black) as a function of the SNR for different val¬ 
ues of the masses. The only case we find where the hybrid 
EoS could be distinguishable from the normal matter one 
is for masses around IAMq. A SNR ~ 30 detection with 
such masses can provide strong BFs in the Jeffreys scale 
if we compare ALF5 to AP4. However, if we compare 
ALF-GCR to GCR we recover smaller BFs. We also find 
that we recover similar results when comparing hybrid 
stars and a normal NS regardless of which one is the cor¬ 
rect star of Nature. This is different from the kaon and 
hyperon cases studied before, where we obtain more con¬ 
clusive BFs when kaons or hyperons are not present in 
NSs. Of course, if the true hybrid EoS of Nature con¬ 
tains weaker strong interactions between the quarks, the 
prospects of detecting a hybrid star reduce even further. 


2. Quark stars 

Our final study case is SQM3 [93], an EoS that contains 
no normal matter at all, but rather it is constructed solely 
with quark matter. Comparing it to a normal matter 
EoS amounts to comparing normal NSs to strange quark 
stars. Although the latter have already been heavily con¬ 
strained [lOlj . it is still interesting to study the bounds 
aLIGO could place on them. However, it is not clear 
what EoS we should compare SQM3 to. Throughout our 
analysis, we compare EoSs that contain exotic matter to 
the EoS that we obtain if we remove the exotic matter, 
but change nothing else in how it is constructed; this is 
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FIG. 6. (Left) m-R relation for the hybrid EoSs and the nuclear matter EoSs they are based on and compared to: AP4 (solid 
black), ALF5 (dashed black), GCR (solid red), and ALF-GGR (dashed red). The pairs of similar symbols denote the injected 
masses. The larger the value of the strong interactions, the larger the deviation between the normal matter and the hybrid 
EoS. (Right) BF for the ALF5/AP4 comparison (black lines) and the ALF-GGR/GCR comparison (red lines) as a function of 
the SNR for different injected masses. BFs are given in favor of the correct model denoted in the legend. When the strong 
interactions between the quarks are close to the value predicted from perturbative calculations, it is possible to distinguish 
between normal and hybrid NSs. 




FIG. 7. (Left) m-R relation for SQM3 and AP4. The pairs of similar symbols denote the injected masses. (Right) BF in favor 
of AP4 (black) and SQM3 (red) as a function of the SNR for different injected masses. BFs are given in favor of the correct 
model. aLIGO will be able to place strong constraints on the existence of strange quark stars, both detecting them if present 
in Nature or strongly disfavoring their existence if not. 


precisely how we defined EoS pairs. In the case of SQM3, 
if we remove the quarks there is no matter left, so there is 
no normal matter counterpart that can form a pair with 
a strange quark star. 

For that reason and in order to arrive at conserva¬ 
tive conclusions, we will compare SQM3 to AP4 [84], a 
reliable EoS that has both not yet been ruled out by 
observations and leads to NSs that are the most simi¬ 


lar to SQM3 quark stars in the set of EoSs we consid¬ 
ered. The left panel of Fig. shows the mass-radius 
relation for SQM3 and AP4, along with symbols that 
indicate injected masses. Being the softest EoS in the 
set of nonruled out EoSs we considered, the AP4 mass- 
radius relation is the farthest to the left in Fig. [T] and 
thus the closest to the SQM3 one. Other normal matter 
EoSs are stiffer than AP4, predicting larger radii for the 
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same mass, which leads to m — R curves even farther 
away from SQM3 than AP4. As we show below, AP4 
is already clearly distinguishable from SQM3, and thus, 
distinguishing between SQM3 and other stiffer normal- 
matter EoSs would be even easier. 

Clearly, SQM3 is very different from all other EoSs 
studied here; its m — R relation differs qualitatively from 
normal matter EoSs, and we thus expect it to be distin¬ 
guishable. Observe that SQM3 fails to produce a 2Mq 
NS, though it is still consistent with the current obser¬ 
vational bound on the maximum NS mass [l()2j within a 
2 -(t statistical error. Of course, it is possible that the star 
observed in |102) is simply a NS, and not a quark star, 
without necessarily ruling out the existence of the latter. 
We have thus decided to study whether aLIGO can dis¬ 
tinguish between strange quark stars and normal NSs. If 
all compact stars that aLIGO sees are NSs or BHs, then 
this would build confidence for the nonexistence of quark 
stars. 

The right panel of Fig. shows the BE in favor of 
AP4 (black) and SQM3 (red) as a function of the SNR 
of the signal for different mass combinations. The really 
high BFs we recover indicate that aLIGO will be able to 
both detect strange quark stars if they exist or provide 
very strong evidence for their nonexistence. For example, 
the detection of a (1.4,1.35)M0 system with SNR = 20 
results in strong BFs in the Jeffreys scale in favor of the 
correct model. For even brighter sources, or less massive 
systems, we recover very strong or even decisive evidence 
in favor of the correct model of Nature, be it strange 
quark stars or normal NSs. 

We have argued that the results given above represent 
the worst case scenario when comparing quarks stars to 
normal NSs, given the reliable normal matter EoSs avail¬ 
able today and not yet ruled out by observations. But 
what if nuclear physicists construct a viable EoS that 
is softer than AP4? Even in this scenario, we can make 
some claims based on the unique shape of SQM3. Revisit 
the left panel of Fig. and notice that between 1.8Mq 
and IMq the radius SQM3 predicts increases by about 
2km. A normal matter EoS that is softer than AP4 will, 
roughly speaking, have a similar shape to AP4, but it will 
be shifted to the left in the m — R plane. Even in this 
case, there will exist some mass between IMq and 1.8Mq 
where the difference between the radius of a quark star 
and a normal NS is at least 1km. Our results indicate 
that typically a radius difference around 0.5km is about 
enough for distinguishability for this type of study, pro¬ 
vided the mass is not close to the maximum mass allowed. 
We can, therefore, claim that even if a softer than AP4 
EoS is constructed, it will still be distinguishable from 
SQM3 for some value of the mass between 1 and I.SMq, 
given an SNR around 30, depending on its exact shape. 



FIG. 8. BF in favor of SV (a kaonless EoS) compared to 
SV222 (a kaon EoS) as a function of mi for m 2 = 1.9Mq 
with SNR= 30. 



FIG. 9. 2D Scatter plot in the mi — m 2 plane for injected 
masses m 2 = 1.9Mq and mi = 1.92M0 and mi = 2.SIM©. 
The red box indicates the injected masses and the vertical 
line gives the maximum mass SV222 can support. 


D. Edge effects 

So far, we have injected signals with masses for which 
both models can produce a NS (or a hybrid star or a 
quark star). In the case of hyperons and kaons, we saw 
that the large value of the masses required to tell EoS 
models apart are close to the maximum mass of the exotic 
matter model. Low strength signals are greatly affected 
by this maximum mass cutoff. Here, we examine the case 
where one of the two objects can only be supported by 
the model that predicts the higher maximum mass. 

We revisit the pair SV-SV222 and fix the SNR to 30 
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and the mass of the smaller NS in the binary to l.QM©. 
We inject signals constructed with SV and in Fig. we 
plot BF in favor of the correct model as a function of 
the NS mass of the larger star mi. The maximum mass 
supported by SV222 is ^ 2 . 11 M 0 . When both simulated 
masses can be supported by the wrong model, then we 
find BFs that are only barely worth mentioning in the 
Jeffreys scale in favor of the correct model. However, 
when mi > 2 .IIM 0 the BF in favor of SV starts increas¬ 
ing and it diverges around mi = 2.35Mq. At this point, 
the chirp mass of the system is so large that it cannot be 
matched by a system with both masses below 2 . 11 M 0 . 
Beyond this point, SV222 cannot produce any systems 
with the correct chirp mass of the injected signal. 

To illustrate this transition we plot the 2—D scatter 
plot of the chain points in the mi — m 2 plane in Fig.j^for 
mi = I. 92 M 0 and mi = 2.34M0. The red boxes mark 
the injected masses and the vertical line is the maximum 
mass SV222 can support. The scatter plots have support 
only along constant chirp mass lines. For mi = 1.92M0 
most points fall on the left of the vertical line, and SV222 
can provide a good match for them. However, as we 
increase mi, more and more points move beyond SV222’s 
maximum mass, resulting in BFs that favor it less and 
less (recall that the BF, as given in Eq. ([^, is the ratio 
of the points in SV222 over the points in SV). When 
mi = 2 . 34 M 0 only a small number of points can be 
supported by SV222, and we recover very strong BFs in 
favor of SV. If we increase the mass even more, no points 
fall on the left of the vertical line, and the BF goes to 
infinity. 


E. Noise curves 

Apart from the Zero-Det., High-P. noise configuration, 
there are a few other tuned noise curves for aLIGO m- 
Among them, the NSNS Optimized configuration gives 
the optimal SNR for a NSNS coalescence, while the 
high frequency one achieves the best sensitivity around 
lOOOHz. Figure [T0| shows these 3 noise curves. The NSNS 
Opt. noise curve is tuned to give the highest SNR by re¬ 
ducing the noise levels around (60 —500)Hz, however this 
comes at the expense of higher noise in the kHz region. 
The finite-size effects that we are looking for manifest 
themselves at frequencies above 400 Hz, suggesting that 
NSNS Opt. is suboptimal for EoS determination. On the 
other hand, the High E. configuration has the lowest noise 
in a small window around lO^Hz, but it is clearly inferior 
for lower frequencies. This loss of SNR at low frequencies 
makes its suitability for EoS studies questionable. 

Eigure [D demonstrates how the standard High-P., 
Zero-Det. configuration is the optimal noise curve for 
measuring the EoS both for systems of constant SNR 
(top panel) and for systems at the same distance (bot¬ 
tom panel). We plot the BF for a system with the same 
parameters (apart from the distance in the top panel) 
achieved with the 3 noise configurations. In the constant 



FIG. 10. Sensitivity curves of various aLIGO configurations. 
The High-Freq. optimized curve (blue dashed line) is very 
sensitive in a small window around 1000 Hz, but it has much 
higher noise at lower frequencies. The NSNS Opt. noise curve 
(red dotted line) has slightly lower noise at frequencies below 
600 Hz, but much higher noise above this. EoS effects become 
important at frequencies above 400 Hz. 


distance case, NSNS Opt. achieves the best SNR value 
while High F. has the worst, as expected. In both cases, 
the lower high-frequency noise of the Zero-Det., High- 
P. configuration together with its low-enough noise at 
lower frequencies outperforms both tuned configurations 
in EoS extraction. This is another manifestation of the 
effect discussed in [5D] ; when we are looking for an effect 
that appears only at certain frequencies, it is not the to¬ 
tal SNR that matters, but the SNR accumulated in those 
frequencies. 


IV. CONCLUSIONS 

We studied whether future GW data from the ad¬ 
vanced, ground-based detectors will allow us to learn 
about the internal composition of NSs. We find that 
aLIGO can efficiently distinguish between NSs with EoSs 
that differ at relatively low to moderate central densities. 
This is the case, for example, for EoSs that model just 
pure quark matter. If the NS is a hybrid of normal and 
quark matter, higher SNR values are required, as well 
as relatively strong interactions between the quarks. On 
the other hand, aLIGO will not be able to efficiently dis¬ 
tinguish between NSs with EoS that differ only at high 
central densities. This is the case, for example, for EoSs 
that model normal matter and those that include a hy- 
peron or a kaon condensate in the inner core. In this 
case, aLIGO would require a very loud detection to be 
able to discern between normal and exotic matter NSs. 

In this paper, we considered kaon, hyperon, and hybrid 
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FIG. 11. BF for the same system with the different aLIGO 
noise curves. In the top panel we keep the SNR constant, 
while at the bottom panel we keep the luminosity distance 
constant. The SNR to which this distance corresponds with 
each noise curve is indicated in the plot. In both cases, 
the Zero-Det., High-P. conhguration gives the highest BFs at 
these distances and SNRs, making it the optimal noise curve 
for EoS studies. 


EoSs with exotic matter parameters within the range al¬ 
lowed from experiments and theoretical calculations, that 
show relatively large difference from the corresponding 
normal matter EoS. However, these exotic matter pa¬ 
rameters can be varied within their experimental range 
to construct exotic matter EoSs that are more similar 
to their corresponding normal matter EoS. If one were 
to marginalize the BFs over such exotic matter parame¬ 
ters by considering as many exotic matter EoSs as possi¬ 
ble with different choices of parameters, calculating the 
BFs and taking the average, the main results in this pa¬ 
per would still hold. In fact, such an analysis would 
strengthen our conclusion since one would find that it 
would be more difficult for aLIGO to distinguish between 
normal and exotic matter EoSs than the results presented 
in the main text. 

Our results could be improved if an accurate descrip¬ 
tion of the merger, where the finite size effects are more 
prominent, is included in the models. For example, the 
effect of hyperons and strange quark matter on GWs 
from merging NS binaries have been studied e.g. in m 
and m respectively. However, NS mergers always oc¬ 
cur at kHZ frequencies]^ where the detector sensitivity 
is lower and it is the inspiral phase that falls in the 
most sensitive frequency band. Moreover, to this day, 
no complete template bank of NSNS mergers exists. It 


would be interesting to study how much information can 
be extracted from the merger phase, and how this could 
aid EoS determination by carrying a full inspiral-merger 
analysis. 

As a side note, we find that the aLIGO noise config¬ 
uration that maximizes the gain of relevant physical in¬ 
formation about NS EoS is the standard Zero-Detuned, 
High-Power one. This is perhaps contrary to the belief 
that a high-frequency tuned noise configuration would 
do best. The reason why the Zero-Detuned, High-Power 
configuration does better is that the finite-size effects 
that depend on the interior composition of the stars ac¬ 
cumulate from roughly a few hundred Hz all the way up 
to merger, and not in a narrow band in the kHz range. 
Of course, these effects are much smaller at hundreds 
of Hz than at contact, however, the noise configuration 
that accumulates the most SNR when finite-size effects 
are non-negligible is the Zero-Detuned, High-Power one. 

The above picture could, again, be altered if the merger 
was included. Numerical relativity studies |l(kimi(i] have 
shown that the EoS leaves an important imprint in the 
merger and post-merger phases in the form of resonant¬ 
like features in the Fourier GW amplitude due to oscilla¬ 
tions of hypermassive NSs. For such features that spike 
in a very narrow band, it may be that a high-frequency 
tuned noise curve would be best. As of today, it is un¬ 
clear how much physical information could be extracted 
from such very high frequency features. 

In our study we attempted to focus on physical ques¬ 
tions of model selection as opposed to analyzing all pos¬ 
sible EoSs constructed. There may exist other interest¬ 
ing physical features of NS interiors, other than the ones 
studied here, that might be worth investigating. To dis¬ 
cover what other features would be interesting to mea¬ 
sure, a stronger synergy between the GW and the nuclear 
physics community should be encouraged. 
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- Appendix A: EoS Classification 

The NS merger frequency cannot be pushed arbitrarily low since In this appendix, we classify various EoSs into Several 

the NS mass cannot exceed roughly 3Mq by causality. categories depending On their physics/compositions and 
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methods used to derive them iminH], and describe 
each class briefly. The basic properties of each class are 
summarized in Table nni 


1. Normal matter EoSs 

First, we focus on EoSs with normal matter, namely 
neutrons, protons, electrons and muons. The EoS can be 
calculated as follows. In the core region, the total energy 
density of the npe^ matter is given by 

e(n„. Up, rie, n^) = eAr(n„, Up) + UnrUnC^ 

+npmpC^ + ee(ne) 

+epinp) , (Al) 

where nt and irii {i = n, p, e, /i) correspond to the num¬ 
ber density and the mass of the npe/r matter respectively, 
ew('«■«: np) represents the total nucleon contribution, cal¬ 
culated from different methods for treating the many- 
body problem with the effective nucleon nucleon interac¬ 
tion (NNI) model constructed based on existing experi¬ 
mental data, ee(ne) and ep^{np) are the energy densities 
of electrons and muons respectively, which one treats as 
free Fermi gases since Coulomb contributions are negli¬ 
gible compared to the kinetic energies. Equilibrium con¬ 
ditions on chemical potentials /i^ (i = n, p,e,/r) for npe/i 
matter with respect to weak interactions are given by 
Mn = Mp + Me 3'iid = /ie with fij = dtjduj. Together 
with charge neutrality, Up = Ue + one can express 
with the baryon density n = + Up. This means that 

one has the total energy density e only in terms of n. The 
pressure can be derived from the first law of thermody¬ 
namics given by 

(la) Variational Chain Summation Method: AP4, GCR 

To calculate the EoS in the nonrelativistic limit, we 
need, in principle, to derive the energy of each component 
of nuclear matter by solving the many-body Schrddinger 
equation. To minimize complexity, Pandharipande [62] 
proposed to estimate the total energy E of the system by 
numerically finding the correlation function in the wave 
function dt that minimizes E = (d>,7J5')/(d>,5'), where 
(,) denotes the inner product and H is the Hamiltonian. 

This variational method with certain constraints on 
the correlation function reduces to solving a simplified 
Schrddinger equation with an interaction potential that 
depends on the baryon density n. Solving this equation 
and performing a cluster series expansion on the Hamilto¬ 
nian to keep up to two-body clusters (namely, neglecting 
three-body clusters and higher), one obtains E{n), which 
is the sum of the kinetic and potential energies. From this 
energy, one can easily derive the pressure using Eq. (|A2|) . 


AP4 [M] uses the Argonne (AV18) NNI model [119] 
and the Urbana IX (UIX) model of TNI [120] . The for¬ 
mer is obtained by fitting the proton-proton and proton- 
neutron scattering data from the Nijmegen group [121] 
and is expected to include all leading many-body corre¬ 
lation effects. AP4 further includes the relativistic IPN 
boost correction to the AV18 model. 

GCR [100] uses the quantum Monte Carlo 
method [122] to systematically study the effect of 
TNIs on the EoS by varying the strength of such inter¬ 
actions and the range of the short-range TNI force. Due 
to the complexity of the method, the EoS is constructed 
for pure neutron matter. The authors use the same 
NNI model as AP4 and effectively take relativistic 
effects into account by assuming that such effects have 
a similar density dependence to TNIs, as shown in [84] . 
Below n < 0.08 fm“^, the EoS is matched to the crust 

EoS [iMns]. 


(Ib) Skyrme-Hartree-Fock (SHF) models: 

SV, SGI, SkI4 

SHF EoSs are constructed by using a Hartree-Fock ap¬ 
proximation with effective, Skyrme type NNI models [63] . 
The ground state energy of nucleon matter is given by 
(4'o|ipAr|'I'o), where is the ground state wave func¬ 
tion that includes nucleon correlations and minimizes the 
energy, while Hat is the bare nuclear Hamiltonian in¬ 
cluding NNIs and TNIs. The mean-field scheme approx¬ 
imates this energy as (<i)o|H|f |4>o) = / 'H{r)d^r, where 
$0 is the Hartree-Fock wave function, is the effective 
nuclear Hamiltonian and H is the Hamiltonian density. 
The standard Skyrme type model uses the parametrized 
Hamiltonian density II26HI28] including both local 
and nonlocal terms, density-dependent terms and a spin- 
orbit interaction. Such a model has up to 10 free pa¬ 
rameters that are determined by fitting to experimental 
data. Apart from these ten parameters, H also depends 
on the neutron and proton local matter, kinetic and spin 
densities, which is given through the Hartree-Fock wave 
functions. 

The SV NNI model [129] uses five out of ten parame¬ 
ters to fit experimental data of the total binding energies 
and charge radii of magic nuclei (forming complete shells 
within nuclei), and it does not include three-body inter¬ 
actions. The SGI NNI model [126] improves previous 
models so that it predicts reasonable values for the in¬ 
compressibility, spin and spin-isospin Landau parameters 
and pairing matrix elements. The SkI4 NNI model [HB] 
includes an additional parameter in the spin-orbit term 
such that it can fit the measurement of charge isotope 
shifts in Sr and Pb isotopes. 


(Ic) Brueckner-Hartree-Fock (BHF) models: 
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Class 

EoS 

Method/Model 

Composition 

Relativistic TNI 

la 

AP4 |84| 

variational 

n, p, e, p 

IPN 

Yes 

la 

GCR dUD] 

variational 

n 

(Yes) 

Yes 

Ib 

SV [85] 

SHF 

n, p, e, p 

No 

No 

Ib 

SGI [85], SkI4 [86] 

SHF 

n, p, e, p 

No 

Yes 

Ic 

DBHF12)(A) [H7j 

BHF 

n, p, e, p 

Yes 

No 

Ic 

MPa [88] 

BHF 

n, p, e, p 

No 

Yes 

Id 

G4 [H2], GA-FSU2.1 [§0] 

RMF 

n, p, e, p 

Yes 

No 


Hb I 

SGI-YBZ6-SAA3 |5T1, SkI4-YBZ6-SAA3 [5T1, 

SHF 

n, p, e, p, H 

No 

Yes 

He 

N1Y5KK* [87], 

BHF 

n, p, e, p, H 

No 

No 

He 

MPaH [88], 

BHF 

n, p, e, p, H 

No 

Yes 

Hd 

H4 [SI] 

RMF 

n, p, e, p, H 

Yes 

No 

HIb 

SGI178 [85] 

SHF 

n, p, e, p, K 

No 

Yes 

HIb 

SV222 [SS] 

SHF 

n, p, e, p, K 

No 

No 

Hid 

GA-FSU2.1-180 [90| 

RMF 

n, p, e, p, K 

Yes 

No 

IVa 

ALF4 [92], ALF5 

variational 

n, p, e, p, TT, Q 

IPN 

Yes 

IVa 

GCR-ALF 

variational 

n, Q 

(Yes) 

Yes 

V 1 

SQM3 [93] 

MIT bag 

Q (u, d, s) 

Yes 

— 


TABLE III. Properties of each EoS class. The upper half corresponds to the EoSs with normal matter only, while the lower 
half corresponds to those with exotic matter. The first column shows the classification number while the second column is an 
example EoSs within each class. The third column describes the method or model that is adopted to calculate each EoS. The 
fourth column shows the composition of matter, where H, K, tt, Q refer to hyperon, kaon, pion and quark respectively. The 
fifth column presents whether relativistic effects have been taken into account (GCR effectively takes relativistic effects into 
account by assuming such effects have a similar density dependence to the short-range TNIs), while the final column shows 
whether three nucleon interactions are included. All the EoSs listed here are valid only at zero temperature. 


DBHF<‘^\A), MPa 


The effective NNI in BHF models is obtained through 
the G-matrix G(n; w) (where uj represents the unper¬ 
turbed energy of the interacting nucleons) by solving the 
Bethe-Goldstone equation self-consistently m- Such 
an equation depends on the bare NNI potential and the 
single particle energy. The only parameters in the the¬ 
ory are those that appear in the former while in the 
nonrelativistic case, the latter is given by the sum of 
the kinetic energy and the single particle potential as 
e(fc) = fc^/2m + U(k\n) with m representing the bare 
nucleon mass. The Brueckner-Hartree-Fock (BHF) ap¬ 
proximation for U{k]n) is given by 

U{k;n)= {kk'\G[n;e{k) + e{k')]\kk')a , 

k'Kkp 

(A3) 


where kp represents the Fermi momentum. The sub¬ 
script “a” refers to the antisymmetrization of the matrix 
element. In the BHF approximation, the energy per nu¬ 
cleon is given by [T30] 


E 

A 


ikl 

10m 


^ ^ U{k-n), 

k<kF 


(A4) 


with A representing the nucleon number. 

In the relativistic model, one introduces a strongly at¬ 
tractive scalar component and a repulsive vector com¬ 
ponent in the nucleon self-energy, which can be self- 
consistently obtained by solving a modified Thompson 


equation [BJ]. Such components add relativistic correc¬ 
tions to the nucleon mass and the single particle energy, 
which then modifies the G matrix and the energy per nu¬ 
cleon. The pressure is again obtained from the first law 
of thermodynamics, given in Eq. (A2). 

DBHF^^)(A) [ST] is constructed by solving the coupled 
Bethe-Saltpeter equations in the rest frame of nuclear 
matter using the Bonn A potential |131] within the “ref¬ 
erence spectrum” approximation, which assumes that the 
values of self-energies are frozen at the Fermi momentum 
at each total baryon density. Nucleons interact through 
exchanges of two scalar (cr, i5), two vector (w, p) and two 
pseudoscalar ( 77 , tt) mesons. 

MPa is constructed within the nonrelativistic 

BHF model with three-nucleon interactions taken into 
account. The NNI is described by the extended soft 
core model m which explicitly includes two-meson and 
meson-pair exchanges. TNIs have two parts, repulsive 
and attractive. In particular, the former is important 
to make the EoS stiffer so that it can support NSs with 
a mass larger than 2Mq. The three-nucleon repulsive 
(TNR) part is modeled by the multi-Pomeron exchange 
potential including the triple and quartic Pomeron ex¬ 
change, where the Pomeron is related to an even num¬ 
ber of gluon exchanges. The strength of the TNR part 
is determined from elastic scattering experi¬ 

ments imiiMj. The three-nucleon attractive part is 


In this paper, we refer to the MPa EoS without hyperons as 
“MPa” and that with hyperons as “MPaH”. 
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added such that it reproduces the nuclear saturation 
property [135] , Having such an interaction model at 
hand, one can calculate the G-matrix, and in turn the 
EoS with chemical equilibrium conditions, charge neu¬ 
trality and baryon number conservation. 


fid) Relativistic mean field (RMF) theory: 

G4, GA-FSU2.1 

Reference [89] constructs the G4 EoS, where the au¬ 
thors model the low-energy strong nuclear interaction as 
the leading exchange of mesons and baryons. One starts 
with the construction of a relativistically covariant La- 
grangian, which includes free baryons, leptons, cr, w and 
p mesons, the tree level meson-baryon interactions and 
perturbative self-interactions for the a meson. The the¬ 
ory is a phenomenological low-energy effective field the¬ 
ory. Such a theory has a small number of parameters, 
which are determined from experiments. The fields are 
replaced by their mean values, assuming that the bulk 
matter is static and homogeneous. The derivatives of the 
fields with respect to time and space also vanish, which 
simplifies the Euler-Lagrange equations. Such equations 
are combined with those of charge and baryon number 
conservations and /3-equilibrium, to yield the Fermi mo¬ 
menta and meson fields as functions of the baryon den¬ 
sity. One then uses these solutions to derive the energy 
and pressure. 

Gupta and Arumugam [^ constructed an EoS us¬ 
ing the effective field theory-motivated RMF (E-RMF) 
model |136] , which has a few additional couplings on top 
of the standard RMF models. The effective Lagrangian 
in the E-RMF model is obtained in a power series of 
fields and their derivatives to a certain order, and all 
the nonrenormalizable couplings are made to be consis¬ 
tent with the symmetries of quantum chromodynamics 
(QGD). The E-RMF model can simultaneously explain 
finite nuclei and infinite matter |137j . In this paper, we 
use the EoS in with the FSU2.1 parameters for the 
nucleon-meson coupling constants, which we call the GA- 
FSU2.1 EoS. 


2. EoSs with exotic matter 

Let us now consider EoSs that also include exotic mat¬ 
ter, like hyperons, kaons and quarks. 


(Ilh) Hyperons with SHF Models: 
SGI-YBZ6-SAA3, SkI4-YBZ6-S 

Reference m included the AA hyperon interactions 
into a few Skyrme-type EoSs, in particular, SGI and 
SkI4. As in constructing SHF models with normal mat¬ 
ter components. Ref. |9Ij parametrized the Hamiltonian 


density for a A hyperon |I38] including the AA inter¬ 
action [139], whose parameters are again determined 
from experiments. Several sets for these parameters ex¬ 
ist [I381II40HI42j . including the YBZ6 model |140j for the 
parameters associated with the NA interaction and the 
SAA3 model [139] for the AA interaction. On top of these 
Skyrme-type interaction potentials, the authors in |91j 
introduced a finite-range force model for the AA interac¬ 
tion. This is because the binding energies calculated from 
the Skyrme-force models are not fully consistent with ex¬ 
perimental data. Parameters in such a non-Skyrme-type 
interaction are determined from the measured binding 
energy of aaHo [143] . 

Based on these interaction potentials, the EoSs are cal¬ 
culated under /3-equilibrium conditions, charge neutrality 
and baryon number conservation. The EoS constructed 
with the SGI (SkI4) NNI model, the YBZ6 model for the 
NA interaction and the SAA3 model for the AA inter¬ 
action is called the SGI-YBZ6-SAA3 (SkI4-YBZ6-SAA3) 
EoS. 


(He) Hyperons with BHF models: 

NIY5KIC, MPaH 

The N1Y5KK* EoS [87] is constructed by extending 
DBHF(^^(A) to include the effect of hyperons (A, 
and S“) with the vector K* and pseudoscalar K mesons 
which induce the baryon-exchange and baryon-transition 
processes. The hyperon-meson coupling constants are 
determined from SU(6) symmetry |I44j . 

The MPaH EoS is an extension of the MPa EoS that 
includes hyperons (A and S“). This EoS is constructed 
under the assumption that TNR-like repulsive interac¬ 
tions among three nucleons work universally for HNN, 
HHN, HHH (as well as NNN), where N and H refer to a 
nucleon and hyperon respectively. 


(Hd) Hyperons with RMF theory: H4 

Reference [HI] constructs EoSs with hyperons using 
RMF theory by adding such effects on top of the G4 
EoS. Hyperons are produced at roughly twice nuclear 
density. The most relevant hyperons are A and hy¬ 
perons, which have the smallest masses. The hyperon- 
meson couplings are taken to be the same among all 
hyperons and are smaller than those between nucleons 
and mesons. The H4 EoS is constructed by setting 
{K,m*/m,Xa) = (300MeV, 0.70, 0.72), where K, m, m* 
and Xu represent the incompressibility, the bare nucleon 
mass, the effective nucleon mass and the ratio between 
the hyperon-meson and nucleon-meson couplings for the 
a meson respectively. 


(HIb) Kaon condensates with SHF models: 
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SV222, SGill8 


For a given density, all the parameters in the nucleon 
and lepton Hamiltonian densities are determined from 
the /3-equilibrium and charge neutrality conditions, and 
the properties of finite nuclei |145) . The kaon Hamil¬ 
tonian density consists of the kinetic and mass terms 
of the kaon and the kaon-nucleon interactions, adopting 
a SU(3) nonlinear chiral model first proposed in |146j . 
One of the parameters in the kaon Hamiltonian den¬ 
sity, which corresponds to the magnitude of the expec¬ 
tation value of the kaon condensate, is determined by 
minimizing the total energy density. The pressure is 
derived from Eq. (A2). The three remaining param¬ 


eters, oi, 02 and 03 in the kaon Hamiltonian density 
are chosen as aims = —67MeV, a 2 ms = 134MeV |147j 
and a^ms = —222MeV for an SV EoS (SV222) and 
asms = — 178MeV for an SGI EoS (SGI178), with 
representing the mass of the strange quark. 


(Illd) Kaon condensates with RMF theory: 

GA-FSU2.1-180 

Reference m also constructed an EoS based on the E- 
RMF theory with kaon (and antikaon) condensates {K~ 
and K'^). On top of the effective Lagrangian for nuclear 
matter, the following Lagrangian for kaon condensates 
are added: Ck = D^K*Df^K—m^K*K, where K = K~ 
or ifO, mK represents the effective mass of kaons that 
depends on the scalar meson fields and is the differ¬ 
ential operator that depends on the vector and isovector 
meson fields. The presence of kaon condensates modi¬ 
fies the expectation values of the meson fields that are 
used to replace the fields themselves in the Lagrangian. 
The total energy density is given by e = e^r + e/c with 
Ck = mKinii- + nj^o), where cat and ck represent the 
energy density of the nucleon phase |148j and kaon con¬ 
densates respectively, while and n^o are the number 
density of K~ and K^. On the other hand, the expres¬ 
sion for pressure is unaffected by the presence of kaon 
condensates [US]. In this paper, we use the EoS with the 
FSU2.1 parameters and the optical potential of a single 
kaon in infinite matter as Uk = — ISOMeV, which we call 
the GA-FSU2.1-180 EoS. 


(IVa) Hybrid nuclear and quark matter 
with variational method: ALF4, ALF5, GGR-ALF 

Reference [S2] constructs the hybrid EoS of the ALF 
family for nuclear and quark matter. The former is mod¬ 
eled using the APR EoS described above. The EoS in the 
low density region is modeled by the standard tabulated 
EoS of [Ml The quark matter EoS is based on 

a physical model, which takes into account both phases 
of normal unpaired quark matter and color-flavor-locked 
(GFL) quark matter. In the latter, quarks form Cooper 


pairs and one-to-one correspondence arises between three 
color pairs and three flavor pairs. The phase with the 
lower free energy is favored. 

Neglecting perturbative QCD corrections, the free en¬ 
ergy (or the grand potential) density of unpaired quark 
matter ilunpio-q) (with qig representing the chemical po¬ 
tential for a quark) contains the kinetic contribution 
from a degenerate free gas of three colors of relativis¬ 
tic quarks and the negative vacuum pressure represented 
by the bag constant B. Such a constant depends on the 
transition density ric at which the quark matter EoS is 
stitched to the nuclear matter one. The authors also in¬ 
troduced a QCD inspired term [72] , which is proportional 
to {1 — c) fig. Here, c is the QCD correction parameter 
where c = 0 corresponds to noninteracting, free quarks. 
Reference |72j carried out a perturbative calculation and 
found c ~ 0.37, although higher order corrections in the 
strong coupling constant may be important for hybrid 
stars, c determines the maximum mass of a star with a 
hybrid EoS and such a mass becomes larger as one in¬ 
creases c. 

Regarding CFL quark matter, the free energy can be 
decomposed into three parts: 

^CFhiRq, Re) = ^OFh’^^iRq) + Me) 

The first term, denotes the contribution 

from the neutral CFL phase. The second term repre¬ 
sents the contribution from pion condensates. Such con¬ 
densates arise when the electron chemical potential ex¬ 
ceeds the mass of the tt~ meson, the lightest negatively 
charged meson in the CFL phase |150H152] . Finally, the 
third term is the contribution from electrons and muons. 
The quark matter EoS is then constructed by using the 
thermodynamic relations: 


p = -n, 


e = — 


d£l 

d In pq 


(A 6 ) 


The matching between the quark and nuclear matter 
EoS is carried out by imposing global charge neutrality 
and pressure balance conditions between the two matter 
phases [15311154j , where the latter condition specifies /ie 
in terms of pq. 

Selected ALF EoSs (ALF 1-4) were used in [118j . In 
particular, ALF4 has (nc,c) = (4.5no,0.3) where ric and 
no represent the transition and saturation densities re¬ 
spectively. The m-R relation for the ALF4 EoS is based 
on the EoS data constructed by the authors in [521 and 
is slightly different from that used in |118j . It seems 
that the latter was constructed by artificially stitching 
the AP4 EoS to the data constructed in [52], but such a 
procedure changes some part of the mixed phase into the 
pure nuclear matter phase, which is energetically more 
disfavored than the mixed phase (so this artificial pure 
nuclear matter phase is unstable). 

Reference m claims that ALF2 can support a 2Mq 
NS (and has a typical NS radius of ^ 13km) but we do not 
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use this EoS due to the following uncertainty. According 
to this EoS has parameters (nc,c) = (3no,0.3) but 
Fig. 6 of [52] shows that the maximum mass with such 
a choice of parametrization is below 2Mq and a typical 
radius is ~ 11km. The ALF2 EoS data used in m has 
a transition of “AP4 —>■ CEL —>■ mixed phase —>■ CEL” 
as one increases the density, although it is not clear if 
such a transition is physically reasonable (ALFl, 3 and 4 
all have a transition of “AP4 —>■ mixed phase —)■ CEL”). 
Instead, following [92], we constructed a new ALE EoS 
(ALF5) with {ric, c) = (2no, 0.4) that can support a 2Mq 


NS. 

We also constructed new hybrid EoSs (GCR-ALF) by 
stitching the nuclear matter GCR EoSs m with the 
quark matter EoSs in |52|- Instead of imposing the 
global charge neutrality condition, we impose charge neu¬ 
trality on each nuclear matter and quark matter (CEL) 
phase |92| . The charge neutrality on the CEL phase leads 
to the absence of electrons (and muons and pion conden¬ 
sates), thus He = 0 and flcFL = 0 = in Eq. (A51. 

The nuclear and quark matter EoSs are matched at a 
chemical potential where the pressure between the two 
phases becomes identical and such a construction leads 
to the absence of the mixed phase [52] ■ As an exam¬ 
ple, Fig. i presents the m-R relation for a GCR-ALF 
EoS with a symmetry energy of Esym = 33.8MeV in the 
nuclear matter EoS (with the crust EoS in |125j l and 
(jic,c) = (2no,0.35) in the quark matter EoS, together 
with the one for the corresponding GCR EoS. We have 
also constructed GCR-ALF EoSs with other parametriza- 
tions, but found that the difference in the m-R relation 
between GCR and GCR-ALF EoSs is typically smaller 
than that shown in Fig [^ 


(V) Strange quark matter: SQM3 


SQM EoSs are constructed from the MIT bag model. 
When one can neglect the mass of quarks, the free energy 
density is given by Q, = CqH ^ -\- B with Cq representing 
a constant. From Eq. (A6|, the EoS for quark matter 


becomes p = 3(e —45). The prefactor “3” changes if one 
includes the effect of the strange quark mass m^. SQM3 
assumes rris = 50MeV and B = 57.39 MeV/fm^, where 
the latter corresponds to the minimum bag constant al¬ 
lowed from the stability of neutrons with respect to a 
spontaneous fusion into strangelets [US]. 


Appendix B: Bayes Factors and Prior Cutoffs 

The NSNS Optimized noise curve has higher noise 
levels at high frequencies than the Zero-Detuned, High- 
Power one. Therefore, as far as EoS determination is con¬ 
cerned, the likelihood calculated with a NSNS Opt. noise 
curve does not contain a lot of information in the high- 
frequency part of the waveform where finite-size effects 
become important. Calculating BFs in a regime where 



FIG. 12. BF in favor of SV compared to SV222 for the system 
with masses (1.95,1.9) with the Zero-Detuned, High-Power 
sensitivity curve (black) and the NSNS Optimized curve (red). 
For high SNR values, where the likelihood dominates, the 
detuning curve gives higher BFs. For lower SNR values with 
the NSNS Opt. curve, we encounter the counterintuitive effect 
of decreasing BFs with increasing SNR values. 



FIG. 13. 2D Scatter plot in the mi — m 2 plane for SNR 10 
(black) and 40 (red) with the NSNS Opt. sensitivity curve 
and with the same system as in Fig. |12| The vertical line 
denotes the maximum SV222 mass. Any points on the right 
of this line necessarily correspond to the SV model. 


it is not the likelihood, but the prior that dominates the 
results can lead to an interesting and rather counterin¬ 
tuitive effect: BFs that initially decrease with increasing 
SNR. In other words, as the signal strength increases, the 
correct model is preferred less and less. 

Figure [^ shows the BF in favor of SV compared to 
SV222 for the (1.95, I.9)Mq system calculated with the 
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Zero-Det., High-P. curve (black) and the NSNS Opt. 
curve (red). The black line is the same as the black 
dashed line of Fig. The red line presents some rather 
interesting behavior. At low SNR, the BF decreases with 
the SNR, while after SNR = 40 it starts increasing, like 
one would expect. In order to understand this effect, con¬ 
sider Fig. where we plot the 2—D scatter plot of the 
chain points in the mi — m 2 plane for the system used in 
Fig. analyzed with NSNS Opt., and SNR 10 (black 
dots) and 40 (red dots). The vertical line indicates the 
maximum mass the competing model SV222 can support. 
This maximum mass is a hard cutoff on the prior mass 
of SV222; no system with masses higher than the max¬ 
imum mass can be produced by SV222, and the model 
reduces to noise. This effectively means that SV222 has 
no posterior weight in that region, and all these points 
correspond to SV. 

As the SNR increases the posterior width decreases. 
At SNR = 40 there are very few points above the mass 
cutoff. That means that with increasing SNR SV222 suf¬ 
fers less and less from this cutoff in its mass prior. This 
explains why the BF in favor of SV decreases as we in¬ 
crease the SNR. Of course, at some values of the SNR 
the differences between the models will start dominating 
over the prior cutoff, and the BF will again start increas¬ 
ing with the SNR. Figure [T^ shows that in this case this 
is true for SNR > 40. 



FIG. 14. BF in favor of hi as a function of the SNR for 
K = 1 (red), K = 0.7 (black), and k = 1.4 (green). The effect 
of the BF decreasing with increasing SNR is present in both 
K = 1 and K = 0.7 cases, but it is more pronounced in the 
case where the two models have a different parameter range. 
As the SNR increases all three lines tend to coincide. This is 
because the likelihood becomes more and more peaked, and at 
some value of the SNR the limits of integration stop affecting 
the integral of the likelihood. 


Appendix C: Toy Model 


To demonstrate the rather counterintuitive effect of 
BFs that favor the wrong model, or decrease with the 
SNR, we construct a simple toy problem, thus obtaining 
a more robust explanation of these effects than the 2-D 
scatter plots presented in the previous appendix. Imag¬ 
ine we receive N data points from a very simple signal 
that obeys d{f) = f and we try to match it with two 
competing 1—D models hi{f) = af and /i 2 (/) = 
where a is the parameter of the models. The likelihood 
for model i is 


N 


L,. = 


\/^c 


:exp<^-^ 


K/)-fe.(/)]^ 

2 ct 2 


(Cl) 


The parameter cr is the standard deviation of the data; 
in GW language, it is 1/SNR. Clearly hi is the correct 
model and Li is maximized when a = 1. On the other 
hand, /12 cannot fit the signal perfectly; L 2 is maximized 
for a = 1.12 with a residual that depends on a. 

Now imagine that the two models have different prior 
ranges for the parameter a. In the case oi hi we have a € 
(0, 2), while /12 allows a £ (0, 2 k ), where k is an arbitrary 
number. The evidence for each model is proportional to 
its likelihood integrated over the range of parameter a, 
while the BF in favor of hi (the correct model) is the 
ratio of the two evidences. Figure shows the BF in 
favor of hi as a function of the SNR (= 1/cr) for three 


different values of k. The case k = 1 corresponds to 2 
models with the same parameter prior range, while k = 
0.7 and k = 1.4 are similar to comparing EoS models 
with different maximum allowed mass. More specifically, 
K = 0.7 corresponds to the case where the wrong model 
has the smaller prior range (the right panels of Figs.|^[^ 
and Fig. 12), while k = 1.4 corresponds to the case where 
it is the correct model that has the smaller prior range 
(the left panels of Figs. and [^. 

The K = 1.4 case is easier to understand. The likeli¬ 
hood of the wrong model ft .2 is integrated over a larger pa¬ 
rameter region a £ (0,2.8) than the correct model where 
a £ (0, 2). When the SNR is low and the models are not 
very different from each other, this can lead to a larger 
evidence for the wrong model. As the signal strength in¬ 
creases the differences between the two models will start 
dominating the evidence and the correct model will end 
up being preferred. Indeed we find that for SNR > 2 the 
BF favors the correct model hi. 

Both K = 0.7 and k = 1 lead to BFs that decrease with 
increasing SNR, however, this effect is more pronounced 
in the k = 0.7 case and it extends to higher values of SNR 
making it easier to identify. When k = 1, one might ex¬ 
pect the Laplace approximation to the evidence to be 
reliable since the posterior width on a is smaller than the 
prior range. However, this is not the case: the likelihood 
for /i 2 is not a Gaussian, and it is this small deviation 
from Gaussianity that we see as a decreasing BF. To vi¬ 
sualize this effect, in Fig. we plot the likelihood for 
hi and /12 for different values of SNR. For low SNR val- 
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FIG. 15. Likelihood for the correct model (solid lines) and 
the wrong model (dotted lines) for different values of the SNR. 
The black vertical line shows the cutoff in the parameter prior 
range when k = 0.7. The likelihood for the wrong model 
does not peak at the injected value of a = 1 and it is not a 
Gaussian. 


ues the likelihoods are essentially the same and BF = 1. 
However, as the SNR increases the area under the red 
dotted line is larger that the area under the red solid 
line, leading to a small decrease in the BF in favor of hi. 
Clearly at sufficiently high SNR the correct model will 
prevail and the BF will start increasing in favor of hi. 


When K = 0.7 the effect of BFs that decrease with SNR 
is stronger and persists for higher values of SNR. Revisit 
Fig. and keep in mind now that the evidence for the 
wrong model is obtained by integrating the likelihood 
up to the black vertical line. As the SNR increases, the 
likelihood for /i 2 is more peaked on the left of the vertical 
line, which means that the area permissible by the prior 
cutoff increases. This leads to an increase of the evidence 
of /i 2 and a BF in favor of hi that decreases. Clearly 
at some point the differences between the models will 
overcome this effect, and the correct model will prevail. 
For our toy model this happens when SNR > 2. 
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